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Abstract 

We study the slow phase of thermally activated magnetic relaxation in finite two-dimensional 
ensembles of dipolar interacting ferromagnetic nanoparticles whose easy axes of magnetization 
are perpendicular to the distribution plane. We develop a new method to numerically simulate 
the magnetic relaxation for the case that the smallest heights of the potential barriers between 
the equilibrium directions of the nanoparticle magnetic moments are much larger than the thermal 
energy. Within this framework, we analyze in detail the role that the correlations of the nanoparticle 
magnetic moments and the finite size of the nanoparticle ensemble play in magnetic relaxation. 
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I. INTRODUCTION 



The role of the dipolar interaction in systems of nanometer-sized ferromagnetic particles, 
or nanoparticle ensembles for short, has been intensively studied in recent years. Such en- 
sembles have numerous technological applications, and it is important to understand their 
magnetic phenomena and processes.^ One of the most complicated problems, where dipolar 
interactions must be taken into account, is thermally activated magnetic relaxation. To de- 
rive the law of magnetic relaxation, i.e., of the dimensionless reduced magnetization, usually 
requires the derivation of the distribution function of the nanoparticle magnetic moments. 
In the simplest case, that of non-interacting nanoparticles with conserved total magnetic 
moments, the distribution function obeys the Fokker-Planck equation.— For nanoparticle 
ensembles with more or less realistic magnetic energy, however, its time-dependent solu- 
tions are not known, and exact results were found mainly for numerical characteristic of the 
relaxation process such as the largest relaxation time . 3 i 4 i 5 i 6 i 7 i 8 

Unfortunately, in the case of dipolar interacting nanoparticles no exact results for the 
magnetic relaxation exist. This fact makes it difficult to check the validity of different ap- 
proximate methods and approaches that are extensively used in this area . 9 ^ ! 11 ^ 2 ' 13 ^ 4 ! 15 ^ 6 ! 17 
The justification of approximations is a very important task because the use of non-rigorous, 
although plausible, approximations can lead to opposite conclusions . 12 > 13 One expects that a 
sufficiently rigorous analysis of the relaxation law can be performed for the simplest systems 
like two-dimensional (2D) ensembles of identical, spherical nanoparticles with conserved 
magnetic moments and large uniaxial perpendicular anisotropy. Such ensembles represent 
an important class of perpendicular magnetic recording media,— and they are convenient 
systems to study experimentally and theoretically the role that the dipolar interaction plays 
in magnetic relaxation. 

Magnetic relaxation in such ensembles was considered first by Lottis, White, and 
Dahlberg 19 within the simplified version of the mean-field approximation. Using the concept 
of a demagnetizing field, the authors wrote down the equation that describes the relaxation 
of magnetization from the initial state, when all nanoparticle magnetic moments are ori- 
ented along a certain direction of the easy axis, to the demagnetized ground state. They 
solved this equation numerically and showed that for a limited time domain relaxation oc- 
curs slower than the Debye model predicts. They approximated the relaxation law p(t) by 
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a stretched-exponential dependence, which, however, does not hold for all times. 

Recently, we studied the influence of the mean and the fluctuating components of the 
dipolar field on the process of magnetic relaxation in those ensembles.— >2i Using the Fokker- 
Planck equation, we derived an equation that describes the so-called slow relaxation, i.e., 
relaxation for times exceeding the time t qe to establish the quasiequilibrium distribution of 
the magnetic moments (t qe ~ 10~ 8 s, see Sec. II B), and we solved it in limiting cases. We 
showed that both the mean and the fluctuating components of the dipolar field enhance 
relaxation, and that for small and large times magnetic relaxation has a Debye character, 
but the corresponding relaxation times can be very different. This difference causes the 
quasi-logarithmic relaxation at intermediate times that was found numerically in Ref. [19 . 

The role that the correlations of directions of the nanoparticle magnetic moments play in 
magnetic relaxation has not yet been clarified. Clearly, correlation effects are very significant, 
and we expect that, due to the antiferromagnetic character of the dipolar interaction in such 
ensembles, they can qualitatively change the relaxation law. The influence of the finite size of 
the nanoparticle ensemble on magnetic relaxation is another important problem, which also 
has not yet been addressed. We expect that, due to the long-range character of the dipolar 
interaction, magnetic relaxation will significantly depend on the ensemble size, especially for 
small times when dipolar fields near the internal and the external magnetic moments can 
be quite different. 

The complexity of these problems forces us to seek numerical solutions. The known 
methods of numerical simulation of magnetic relaxation, such as directly integrating the 
stochastic Landau-Lifshitz equation , 22 i 23 i 24 the conventional Monte Carlo method , 25l2e and 
the time quantified Monte Carlo method 2 ^ are not suitable for our purposes. The main rea- 
sons are the following. To integrate the Landau-Lifshitz equation, the integration time step 
must be smaller than the inverse of the precession frequency of the nanoparticle magnetic 
moments (~ 10~ n s). Therefore, this method usually works only for the description of the 
fast magnetic relaxation, i.e., relaxation on time scales smaller than t qe . The conventional 
Monte Carlo method is not suitable, since each Monte Carlo step has no physical time as- 
sociated with it. The time quantified Monte Carlo method also cannot be applied to our 
situation; the number of Monte Carlo steps that are necessary to calculate the relaxation 
law on times comparable with the relaxation time becomes prohibitively large in the case 
of high potential barriers between the equilibrium directions of the magnetic moments (see 
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Sec. Ill B). Further, that method is valid only in the high damping limit, i.e., if there is no 
precession of the magnetic moments. 

In this paper we develop a new method to numerically simulate thermally activated 
magnetic relaxation in finite 2D ensembles of dipolar interacting ferromagnetic nanoparticles. 
We consider the case where the nanoparticles with uniaxial anisotropy occupy the sites of 
a square lattice and their easy axes of magnetization are perpendicular to the lattice plane. 
We develop an equation that relates the magnetization of the ensemble at the next time 
step to the known state of the nanoparticle ensemble at the previous time step, and a 
numerical procedure that defines the ensemble state at the next time step. To derive the 
probability densities for the reorientation of the nanoparticle magnetic moments, contained 
in this equation, we exploit that they can be represented via the mean times for magnetic 
moments to reorient, or, in other words, via the so-called mean first-passage times, and 
calculate these times using the backward Fokker-Planck equation. 

The paper is organized as follows. In Sec. II, we introduce the equation mentioned above 
and derive rigorous expressions for the probability densities of reorientation of the nanopar- 
ticle magnetic moments. The algorithm for the numerical calculation of the relaxation law 
is described in Sec. III. In the same section we present the numerical results and analyze the 
features of the magnetic relaxation caused by both the correlations of the nanoparticle mag- 
netic moments and the finiteness of the nanoparticle ensemble. We summarize our results 
in Sec. IV. 

II. ANALYTICAL RESULTS 

We consider a system of N uniaxial and identical spherical ferromagnetic nanoparticles 
with a radius r. We assume that the nanoparticle centers occupy the sites of a square lattice 
of size Ld x Ld \{L + l) 2 = N] and lattice spacing d(> 2r). The easy axes of nanoparticles 
magnetization are perpendicular to the lattice plane (xy-plane), and at the initial time t = 
all magnetic moments mi(t) (the index i labels the nanoparticles) are oriented along the z- 
axis (see Fig. 1). We also assume that the smallest heights AUi of the potential barriers 
between the equilibrium directions of the nanoparticle magnetic moments are much larger 
than the thermal energy ksT {ks is the Boltzmann constant, T is the absolute temperature), 
i.e., the condition Si = AUi/ksT 3> 1 holds for all nanoparticles. The main goal of this 
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section is to find the relation between the reduced magnetization at times t and t + r. 
A. Equation for the reduced magnetization 

For £j ^> 1, the vectors nij(t) fluctuate within small vicinities of the positive and negative 
directions of the z-axis, and they are reoriented only rarely. Consequently, the average 
numbers of positively and negatively oriented magnetic moments have well-defined values 
N + (t) and iV_(i), respectively, at any instant t. Since the number of magnetic moments that 
at time t have reoriented is much less than N, the approximate relation N + (t) + N_(t) ~ iV 
holds, and we can define the reduced magnetization of the nanoparticle ensemble as p(t) = 
2N + (t)/N — 1. Let us define also the state of that ensemble. We assume that the state 
of the nanoparticle ensemble at time t is known if the directions of all magnetic moments 
are known, i.e., we describe the ensemble state by the set of signs <Ti(t) = &i (i — 1, ■■■,N), 
where <j{ = + or — depending on whether the vector nij(t) fluctuates around the positive or 
negative direction of the z-axis. 

Given the ensemble state, neglecting the fluctuations of nij(t), and taking into account 
that approximately rrij(t) = <Jime z for the time intervals between the reorientations, we can 
write the local dipolar field hj(t) acting on the magnetic moment m,(t) as hj(t) = hi(t)e z . 
Here 



m = \m.i(t)\, e z is the unit vector along the z-axis, and is the distance between the 
centers of corresponding nanoparticles. If at time t the magnetic moments do not undergo 
reorientations, then each nanoparticle is under the influence of the local dipolar field (j2.1j) . 
Even if some magnetic moments are reoriented, their number is much less than N because 
Si 1, and formula (j2.1|) remains approximately valid. For sufficiently small times intervals 
we can consider therefore the ensemble of interacting nanoparticles as a system of indepen- 
dent magnetic moments, each of which feels its own external magnetic field h$(t). This fact 
significantly simplifies the numerical investigation of the magnetic relaxation in ensembles 
of dipolar interacting nanoparticles. 

Let us assume that the probabilities of reorientation per unit time w a (t;j) (i.e., the 
probability densities of reorientation) of the vectors n%(t) (J = 1, ...,N) from the positive 
direction of the z-axis (if <7j = +) and from the negative one (if <jj = — ) are known. We 
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also assume that on the interval (t, t + r) the probabilities of two and more reorientations 
of mj(t) are negligibly small. Then, taking into account that N + (t + r) — N + (t) is equal to 
the difference between the number of reorientations from the negative direction of the z-axis 
and the number of reorientations from the positive direction of the z-axis, we obtain 

2r N 

p(t + r) - pit) = -— E ^.,<'-r, ,[!:.,) . (2.2) 

iV 3=1 

The probability densities w aj (t;j) depend on the local field hj(t), and Eq. ()2.2j) can be 
applied if the ensemble state at time t is known. However Eq. (|2.2j) is not an iterative 
equation for the ensemble state; it only defines p(t + r) but not the ensemble state at time 
t + r. In order to use Eq. ()2.2j) as the recurrence equation for the calculation of the law of 
magnetic relaxation, we need to determine the values w aj (t;j) and develop a procedure to 
find the state of the nanoparticle ensemble at time t + r, if its state at time t is known. We 
will describe that procedure in the next section. Below we calculate the probability densities 



B. Probability densities of reorientation 

The probability densities of reorientation are given by w a ,(t;j) = 1 / 'tV '(£; j) , where 
t s 3 (t; j) are the mean times that the magnetic moment nij(t) spends pointing in the positive 
(when Oj = +) and the negative (when aj = — ) directions of the z-axis. These times can be 
represented as tV(t\j) = 2tm(t;j), where tm(t;j) are the mean times for Tsij(t) to reach for 
the first time the state with a maximum value of the nanoparticle magnetic energy Wj. The 
factor 2 takes into account the fact that from that state the magnetic moment nij(t) can 
transit to the state Oj = + or <jj = — with probability 1/2. In our case, the magnetic energy 
Wj includes the anisotropy energy —(H a /2m)mj Z (t) and the Zeeman energy —hj(t)mj Z (t), 
so that it has axial symmetry and 

Wj EE Wj(9j(t),t) 

= -^H a m [cos 2 6j(t) + 2bj(t) cos^-(t)] . (2.3) 

Here H a is the anisotropy field, 9j(t) is the polar angle of mj(t), and bj(t) = hj(t)/H a 
(assuming that two equilibrium directions exist for each magnetic moment, < 1 for all 
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nanoparticles) . Accordingly, the state corresponding to the maximum value of Wj is defined 
by the polar angle 

Qj(t) = arccos [-bj(t)] . (2.4) 

From the mathematical point of view, the calculation of the mean times tm(t;j) is a 
particular case of a general problem, known in the theory of Markovian processes as the 
first-passage time problem.^ 8 * In our case, the Markovian process is the vector m.j(t), and 
the level set of first passages for m.j(t) is the conical surface defined by Eq. (|2.4)1 . We describe 
the dynamics of the nanoparticle magnetic moments xtij(t) = m,, by the system of stochastic 
Landau-Lifshitz equations 

thj = -7m.j x (H 3 - + iij) - ^ x (m 3 x H,), (2.5) 

fit 

where j = 1, N, j(> 0) is the gyromagnetic ratio, A(-C 1) is the damping parameter, 

H, = = H a [cos9j(t) + b 3 (t)]e z (2.6) 

is the effective magnetic field acting on uij, and n, = xij(t) is the thermal magnetic field 
that models the action of the thermostat. The thermal field is approximated by Gaussian 



white noise with zero mean values tij(t) = [the overbar denotes averaging with respect to 
the sample paths of tij(t)] and correlations functions 



ni a (ti)njp(t 2 ) = 2A5 i:j 5 a p5(t2 - (2.7) 

Here n ia (t) (a = x,y,z) are the Cartesian components of iij(t), A = XksT/'ym is the 
intensity of the thermal magnetic field, 5ij is the Kronecker symbol, and 5(t) is the Dirac 5 
function. 

If we treat the local dipolar fields hj(t) as external magnetic fields, then we can consider 
the nanoparticles to be independent. In other words, in this case the stochastic Landau- 
Lifshitz equations (|2.5J) are independent, and the dynamics of each magnetic moment is 
described separately. Let Pj = Pj($j,t\$j,t') be the conditional probability density that 
9j(t) = $j given that 9j(t') = d'j (t > t'). [Note that in the case of axial symmetry Pj does 
not depend on the azimuthal angle of m^.] Then, using the Stratonovich interpretation^ of 
Eq. (|2.5|) and applying standard methods, 28 we can write for Pj the forward Fokker-Planck 
equation 



dPj d 



dt d$j 



A 7 dWji&jtt) 



m dftj 



- AY cot tfj 



P., 



,d 2 R 

n 

and the backward Fokker-Planck equation 



+A7 W (2 ' 8) 



dt> 



d 2 P- 



A7 2 cot 



~ AY W (2 ' 9) 

As a rule, the study of the magnetic properties of nanoparticle ensembles is based on 
forward Fokker-Planck equations similar to Eq. ()2.8j) . which allow us to express the sta- 
tistical characteristics of ensembles as functions of time t. At the same time, backward 
Fokker-Planck equations are very useful to describe the thermally induced reversal of the 
nanoparticle magnetic moments." We use the backward Fokker-Planck equation ()2.9j) to 
calculate the mean first-passage times tm(t;j). 

To use Eq. ()2.2j) as the recurrence equation for finding the reduced magnetization at the 
discrete times t — t n (n = 0, 1, M, t = 0, t n+ \ > t n ), we need to calculate tm (t n ; j) for n = 
0, 1, M — 1. Since to each time t n corresponds the angle Qj(t n ), it is necessary in Eq. ()2.9|) 
to replace Wj^pt') by W^(^-,t n ). In other words, to find tm(t n ;j) we must use Eq. ()2.9|) 
with an energy term Wj($j,t') that does not depend on t' . This important requirement 
results in a condition of homogeneity for the random process 9j(t), Pj{i9j, t\$'j, t') = Pji^&j, t — 
t'l'dpO) and significantly simplifies the problem. 

To calculate tm(t n ;j), we first introduce the mean times Tj = Tp("dpt n ), the time 
necessary for Bj(t) [6j(0) = e (0,Qj(t n )) if Oj = +, and ^ G (Qj(t n ),ir) if tr, = -] to 

first reach the angle Qj(t n ). The desired times are expressed through Tj as 

Q(t n ;j)=Tp(n(l-a 3 l)/2;t n ), (2.10) 

and the values Tj themselves are represented in the form 

roo /•7r(l-cr J l)/2+%(t n )(l+crjl)/2 

T = du dd 

JO Jn j (t„)(l-cr J l)/2 

xPj(d,u\^,0). (2.11) 

Taking into account the initial condition Pj($j, 0\$'j, 0) = — the homogeneity condi- 
tion Pj($j, t\"d'j, t') = Pj($j,t — t'l'd'j, 0), and the expression ()2.3|) . we obtain after integration 
of both sides of the modified equation ()2.9j) over u — t — t' and 1? = ■dj as in Eq. (j2.11|) the 
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ordinary differential equation for Tj 

d 2 T- dT- 

+ [cot 7?; - 2a(6 i (0 + cos^) sm#J]^ = -at r (2.12) 

(a = H a m/2k B T, t r = 2/X 1 H a ). 

To find the unique solution of Eq. (|2.12j) . we need to impose two boundary conditions 
for the mean times Tj. The first condition follows immediately from the definition of these 
times: Tj\$i. = Q.(t n ) = 0. We can find the second by analysing the solutions of Eq. (j2.12|) for 
small vicinities of the angles ^ = and i?^ = n. There Eq. (|2.12|) is reduced to 

d 2 Tj 1 dTj 

Wf + Wp tt(1 - <x,l)/2 W 3 ~ ~ atr ' (2,13) 
and its general solution is given by 

Tj = Cj In - tt(1 - (Tj\)/2\ + dj 

-at r [^.-7r(l-a,l)/2] 2 /4, (2.14) 

where Cj and dj are constants of integration. Since Tj are bounded quantities, the condi- 
tion Cj = must hold. This condition can be represented equivalently in the form of the 
second boundary condition: cfZ^/d^l^^i-^ij/a = 0. Note that these boundary conditions 
correspond to the case where so-called absorbing and reflecting barriers^ are placed at the 
points i?'. = ftj(t n ) and fi'j = 0, ir, respectively. 

Solving Eq. (|2.12j) with these boundary conditions by the method of variation of 
constants^ and using the representation (|2.1()jl . we obtain the rigorous formula 

„1 e -a[x+a-jbj{t n )] 2 



/I g— up 
dx 



&jbj(t n ) 1 — X 2 

1 dye a[y+a ^ {tn)] \ (2.15) 

which is valid for arbitrary a and \bj{t n )\ < 1. Using Eq. (|2.15J) . let us calculate w a ,(t n ;j) 
for £j > 1. According to Eq. Ij'i.Hj) . the heights A£/j- 3 (t) of the potential barrier between the 
equilibrium directions of m., can be written in the form AUj 3 (t) = ~H a m[l + <jjbj(t)] 2 , and, 
since AUj = mmAUj 3 (t), the condition Sj ^> 1 leads to a[l + ajbj{t n )] 2 ^> 1. Taking into 
account that the asymptotic formulas 



hold as a[l + crjbj(t n )} 2 — > oo, we find in the same limit 

^(tn;j) = ^y|[l-6j(t„)][l+ ( 7 i 6 i (t ri )] 

xe -«[i+ ff3 t 3 (t„)p ^ 216 ^ 

Note that Eq. ()2.16|) follows also from Brown's results^ for isolated nanoparticles in a lon- 
gitudinal external field obtained with the forward Fokker-Planck equation for Ej ^> 1. We 
have presented here an alternative derivation of Eq. (j2.16|) based on the backward Fokker- 
Planck equation, because within this approach the mean first-passage times tm (t; j) and the 
probability densities of reorientation w a .(t n ]j) are calculated exactly for arbitrary Ej. 

If the condition Ej ^> 1 holds for all nanoparticles and the ensemble state at time t = t n 
is known, then for the same time we can find the dipolar fields acting on each nanoparticle, 
using the formula (|2.1jl . and calculate the probability densities of reorientation of each 
magnetic moment, using the formula ([2. 16)1 . 



C. Mean-field approximation 

To illustrate the influence of the correlations of the magnetic moments and of the finite 
size of the nanoparticle ensemble on the magnetic relaxation, we must first calculate the 
relaxation law p m f(t) for an infinite lattice within the mean-field approximation. To this 
end, we derive the equation that this relaxation law satisfies, based on the results obtained 
above. Since within the mean-field approximation the same mean dipolar field acts on all 
magnetic moments, it is necessary in Eq. (J2.16)) to replace bj(t) (we drop the index n in t n ) 
by b(t) = bj(t). This implies that all magnetic moments for which <jj = + and all magnetic 
moments for which <jj = — are reoriented with the same probability densities, w + (t) and 
W-(t) respectively, where 

w ± (t) = ^[1 - 6 2 (t)][l ± b(t)] e -<*±W (2.17) 

t r V 7T 

The function b(t) is given by^ 

in 

b(t) = -9-034— p mf (t), (2.18) 

therefore the probability densities w±(t) depend on t only via the reduced magnetization 
p m f(t), i.e., w±(t) = w±(p m f(t)). Finally, using the equality J2j &j = N + (t) — N_(t) and the 
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definition of p(t), we obtain from Eq. ()2.2|) for r — > and N ^ oo the required differential 
equation 

Pmf (t) = -Pmf (t) [W+ (t) + W_ (*)] ~ W + (t) +W-(t) (2.19) 

(Pm/(0) = I)? which defines the law of magnetic relaxation in the mean-field approximation. 
Note that the same equation follows from the solution of the forward Fokker-Planck equation 

Calling the right hand side of Eq. (|2.19|) —F(p m f(t)), we can reduce this equation to the 
integral form 

' iV t. (2.20) 



l Pmf (t) F(x) 

Its solution for small and large times yields 2 - Pmf if) = 1 — t/r and p m f(t) oc exp(— t/Voo), 
respectively, where 

r = t r J ; — (2.21) 

V a 4(1 -£ 2 )(1 -£) V ; 

is the initial relaxation time, 

fW e a 

Too = t r \ ; ; r^r (2.22) 

V a 4[1 + (2a - 1)^] V ; 

is the final relaxation time, and £ = —6(0) (0 < £ < 1) is a parameter characterizing the 
intensity of dipolar interaction on an infinite lattice. According to Eqs. (|2.21j) and (|2.22j) . the 
relaxation process in ensembles of dipolar interacting nanoparticles is approximately char- 
acterized by two relaxation times, To and r m , while in the case of non-interacting nanoparti- 



cles, i.e., £ = 0, it is characterized by the single relaxation time r n = t r y ir/16aexp a. Since 
r n > t and r n > r^, the dipolar interaction enhances relaxation, and since r < r^, the 
relaxation rate decreases with time. For ensembles where the value of £ is not too small, the 
strong inequality r <^ usually holds, and the decrease can be very large. 

Note that the description of magnetic relaxation based on Eqs. (J2.2j) and (J2.19j) is valid if 
the quasiequilibrium distribution of the nanoparticle magnetic moments is established, i.e., 
if2£ t > t qe ~ at r . In other words, these equations describe the slow phase of magnetic re- 
laxation. For t ~ tq e , the probability of reorientation of the nanoparticle magnetic moments 
from the initial state is vanishingly small. Therefore we can transfer the origin of time to 
an arbitrary point t ~ t qe and, since for a » 1 and t ~ t qe the approximate equalities 



rrijz(t) ~ m hold, use the initial conditions p(0) = 1 and p m /(0) = 1. 
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III. NUMERICAL SIMULATIONS 



A. The computational algorithm 

According to the results of the previous section, to compute the law of magnetic relaxation 
in some time interval (0, t^) it is necessary to know the states of the nanoparticle ensemble 
at the discrete times t = t n (n = 0, 1, ...,M — 1). The state for n = 0, i.e., for t = 0, is 
known from the initial condition: o"j(0) = + for all j. To find the state at any other time 
we proceed as follows. First we assume that the state of the nanoparticle ensemble at the 
time t = t n is known. This means that the set A + (t n ) of numbers j for which aj(t n ) = +, 
and the set A_(t n ) of numbers j for which <Jj{t n ) = — are fully defined. It is evident that 
the set A + (t n ) contains N + (t n ) elements, and the set A_(t n ) contains N_(t n ) elements. 

Next, assuming that the time interval At n+1 = t n+ i — t n is small enough, we introduce 
the average numbers of reorientations 

v±{t n ,t n+1 ) = At n+1 w ±( t n,j), (3.1) 

jeA±(t n ) 

that occur during At n+ i for the sets of positively (upper sign) and negatively (lower sign) 
oriented magnetic moments. Strictly speaking, Eq. (|3.1|) is valid if the strong inequality 
At n+ i ma,x{w±(t n ; j)} <^ 1 holds. Its use can drastically increase the time required for the 
computation of the relaxation law in some cases. Therefore, instead of the exact represen- 
tation (jH.lj) we use the approximate one 

v±(t n ,t n+ i) = U(At n+ iW±(t n ;j)) (3.2) 

ieA±(t n ) 

[U(x) = x if x < 1, and U(x) = 1 if x > 1], which is valid if the weaker condition 
v±(t n ,t n+ i) <^ N holds, and from Eq. ()2.2J) we obtain 

2 

P(t n +l) = P(t n ) ~ - Kfc.Wl) - f-(t n> *n+l)]- (3-3) 

Equations (|3.3|) . ([3.2)1 . (|2.16|) . and (|2.1|) allow us to calculate the reduced magnetization 
at time t = t n+ i, if the nanoparticle state at time t = t n is known. To find the nanoparticle 
state at time t = t„+i, we need to choose the sites where the magnetic moments must be 
reoriented. To reflect the random character of the thermal fluctuations, these sites should be 
chosen randomly, while at the same time preference should be given to those sites that have 
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larger probabilities of reorientation. To satisfy both requirements we proceed in the following 
way. First we choose the time steps At n+ \. Since the number of magnetic moments that 
are reoriented per unit time can appreciably decrease with time, we select steps of varying 
length, At n+ i = r)[w + (t n ) + The parameter r] must be chosen small enough to 

satisfy the condition u±(t n , t n+ i) <C iV (in our calculations r\ = 5 x 1CT 3 ). Then we calculate 
the values At n+ iw + (t n ; j) for j G A + (t n ), and using the formula ()3.2|) we find the average 
number of reorientations 

i) =r+(t n ,t n+ i) + At n+1 w + (t n ;j), (3.4) 

that occur during the time interval At n+1 in the set of positively oriented magnetic moments. 
Here r + (t n ,t n+ i) is the number of lattice sites where At n+ iw + (t n ; j) > 1, and A' + {t n ) is 
the set of lattice sites where At n+ iw + {t n ; j) < 1. Further, we introduce the number of 
reorientations as n + (t n ,t n+ i) = [u + (t n , t n +i)] + 1, where [v + (t n , i n +i)] is the integer part of 
v + (t n ,t n+ i), and / = or 1 with probability p = u + (t n , t n+1 ) - [v+(t n , t n+ x)} andpi = l-p , 
respectively. Using a random number generator, we obtain a value for n + (t n , t n +i)- 

Among the n + (t n ,t n+ i) magnetic moments that must be inverted at time t = t n+ i, we 
immediately invert the r + (t n ,t n+ i) magnetic moments at lattice sites where the condition 
At n+ iw + (t n ; j) > 1 holds. [Recall that a one-to-one correspondence exists between the 
lattice sites and numbers j.} To find the remaining n + (t n , t n+ i) — r + (t n ,t n+1 ) lattice sites 
where the magnetic moments have to be inverted, we first generate a random number that 
lies in the interval of length J2jeA' + (t n ) w +(tn', j)- This interval contains N + (t n ) — r + {t n , t n+ i) 
subintervals of lengths w + (t n ; j). We store the number j of the subinterval (i.e., the position 
of the site) that contains the random number in memory, and then that subinterval is 
removed. Next we generate a random number that lies in the new interval formed by the 
remaining subintervals. The number j of the subinterval that contains this random number 
is again stored in memory, and then this subinterval is also removed. Iterating this procedure 
n + {t n , t n+ i) — r + (t n , t n+ i) times, we find all n + (t n , t n+ i) lattice sites where positively oriented 
magnetic moments must be inverted at time t = t n +i- 

Introducing in the same way the average number of reorientations z/„(t n , t n+ i) that occur 
in the set of negatively oriented magnetic moments, and using the procedure described above, 
we determine n_(t n ,t n+1 ) lattice sites where these magnetic moments must be inverted at 
time t = t n+ \. Since the ensemble state at time t = t n is known, the ensemble state at 
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t = t n+ i, i.e., after the inversion of n+(t n ,t n +i) + n_(t n ,t n+ i) magnetic moments on well 
defined lattice sites, is known too. Taking the latter state as the initial state, we can find in 
the same manner the ensemble state at time t = t„+2, and so on. 

Using the known state of the nanoparticle ensemble at the initial time t — and applying 
the algorithm described above, we can find the states for all times t = t„ (n = 1, M — 1). 
Since our algorithm is a probabilistic one, the reduced magnetization calculated by the for- 
mula ()3.3|) is a random quantity. Let us designate that random reduced magnetization in the 
kth numerical experiment as p k sim {t n ). [A numerical experiment consists of one application 
of the algorithm to determine the ensemble states at all times t = t n .} Then we define the 
numerically simulated relaxation law as 



where K is the number of numerical experiments. To avoid any misunderstanding, we 
emphasize that within the proposed algorithm the dipolar field (|2.1|) is calculated exactly, 
and it is re-calculated after each time step. 

B. Numerical results and discussion 

We have used our analytical results and the numerical algorithm described above to 
study the role that the finite size of the nanoparticle ensemble and the correlations of the 
nanoparticle magnetic moments play in magnetic relaxation. We found that the reduced 
magnetization p S i m (t) (t > 0) decreases, when the parameter L, a measure of the ensemble 
size, increases, i.e., p S im(X)|zi > Psim(f)\L 2 if ^2 > Li, and p S i m (t)\L tends to the limiting value 
Piimif) as L — > oo. We explain such behavior of p S im(t) as follows. Increasing L leads to an 
increase, on average, of the local dipolar fields acting on the nanoparticle magnetic moments. 
As a result, the average of the probability densities of reorientation of the positively oriented 
magnetic moments increases, and the average of the probability densities of reorientation of 
the negatively oriented magnetic moments decreases. According to Eq. (|2.2|) . this means 
that p S im{t) decreases when L grows. 

To verify this statement, we have calculated p S i m (t) for different ensembles of Co nanopar- 
ticles characterized by the parameters H a = 6400 Oe, m/V = 1400 G (V is the nanoparticle 
volume), A = 0.2, and r = 4 nm. As an illustration, the function p s i m (t), obtained at 




k=l 



(3.5) 
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T = 300 K, d = 3r, L = 50, and K = 100, and the approximate function pii m {t) are shown 
in Fig. 2. We found the latter function in the same way as p S im(t), but, to exclude boundary 
effects, we assume that the basic nanoparticle ensemble (for which we chose L = 100) is 
surrounded by eight identical ensembles, and each nanoparticle from the basic ensemble is 
considered as a central one in the square box of the same size (i.e., L = 100) and interacts 
only with the nanoparticles which belong to this box. In Fig. 2, we also show the func- 
tion p m f{t) calculated via the numerical solution of Eq. ()2.19|) for an infinite ensemble of 
Co nanoparticles with the same parameters. Note that in this case a ~ 29.01, £ ~ 0.31, 
t r « 8.85 x 10^ n s, r « 1.33 x 10~ 5 s, « 1.56 s, and r n « 28.89 s. 

Since at t = the local dipolar field for an infinite ensemble is always larger than the 
highest local dipolar field for a finite one, the condition p s i m (t) > p m f(t) (t > 0) must hold 
for small enough times. We expect that the same condition holds also for large enough times, 
since correlations of the nanoparticle magnetic moments lead to slower magnetic relaxation 
in the final phase than the mean-field theory predicts. As to the relation between p S im{t) 
and p m f(t) at the intermediate times, its character at a fixed temperature depends on the 
ensemble size, i.e., on the parameter L. 

To explain this dependence, we note first that at small times magnetic relaxation for 
finite nanoparticle ensembles occurs faster than in the case where the local dipolar fields are 
replaced by their average value, i.e., the mean-field approximation. Indeed, in the initial 
phase of magnetic relaxation only a small number of the nanoparticle magnetic moments 
is subjected to reorientation. In this case, the reoriented and most of the non-reoriented 
magnetic moments are under the action of the local dipolar fields, which exceed the mean 
dipolar field. This means that w_(t;j) < w_(t) for j E A_(t), w + (t;j) > w + (t) for most 
j G A + (t), and therefore the actual magnetic relaxation occurs faster than the mean- field 
approximation predicts. [We emphasize that this conclusion is valid for the initial phase 
of magnetic relaxation for finite as well as infinite nanoparticle ensembles.] Furthermore, 
taking into account that an increase in the size of the nanoparticle ensemble leads to an 
increase, on average, of the local dipolar fields, we expect the following behavior for the 
dependence of p S i m {t) on L (for an illustration, see Fig. 2). If in the nanoparticle ensemble 
the highest local dipolar field at t = is small enough in comparison to the case of an 
infinite ensemble, i.e., if the parameter L does not exceed the critical value L cr = L cr (T), 
then psimit) > p m f(t) for all t > (curve 1 in Fig. 2). At L — L cr the curves p S i m (t) and 
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p m f(t) have a tangency point, and for L > L cr they intersect at times t = tu n and t = tn n 
(curve 4, tu n ~ 2.25 x 10~ 5 s, t^m ~ 2.09 x 10 -2 s). As L is increased, the time ti in of the 
first intersection decreases, and the time t 2 i n of the second one increases. As a result, for 
L — ► oo we have p S im{t) —>■ pu m (t), tu n —>■ 0, and t 2in tends to the limiting value t in (curve 
2, t in ~ 0.46 s). 

To characterize the difference between p S i m (t) and p m f(t), we introduce the parameter 
Xl{P) — [Psim{t)—pmf(t)]/psim(t)- Its dependence on t for the same ensembles of Co nanopar- 
ticles is shown in Fig. 3. The nonzero value of XzXO is caused by both the finite size of the 
nanoparticle ensemble and the correlations of the nanoparticle magnetic moments. Correla- 
tions significantly change the relaxation law, and their role grows with time, i.e., Xi{t) — ► 1 
as t — > oo. 

The fact that the probability densities of reorientation w a (t;j), Eq. ()2.16|) . depend expo- 
nentially on the large parameter a has two consequences. The first is obvious, namely, the 
relaxation law p S i m (t) strongly depends on temperature due to the inverse proportionality 
of a on T. The second is more complicated and refers to the time dependence of p S i m (t) 
and Pmf(t) for different T. According to the previous results, if at a certain temperature 
the parameter L satisfies the condition L < L cr , then p S i m (t) > p m f(t) for all t > 0. As T 
decreases, the probability densities w aj (t;j) decrease with different rates, and the smaller 
the temperature becomes, the more their relative values differ. This means that as T is 
reduced, the reorientation of the nanoparticle magnetic moments predominantly occurs at 
sites where w aj (t;j) are the largest. As a consequence, for small times the difference be- 
tween p S im{t) and the relaxation law derived by the mean- field approximation grows as T 
decreases. Therefore, if at a given temperature the condition L < L cr holds and the values 
of L and L cr do not differ too much, then the curves p S im(f) and p m f(t) can intersect at 
smaller temperatures. The plots of p S i m (t) calculated for ensembles of Co nanoparticles for 
L = 50 and T = 300 K (see Fig. 2), and for L = 50 and T = 150 K (see Fig. 4) demonstrate 
this statement. In the latter case calculations yield a w 58.02, tq w 10.72 s, r M w 2.24 x 10 12 
s, r n w 8.11 x 10 13 s, t Un ~ 56.12 s, and t 2in ~ 6.83 x 10 11 s. 

The relaxation laws calculated above can not be determined using the Monte Carlo 
method with time step quantification. According to Ref. |27|, the time interval At that 
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corresponds to one Monte Carlo step is written in our notations as 



At 



R 2 (l + X 2 )m 
20k B T\~f 



(3.6) 



(R < 1), and the number M = fj,r n /At of the Monte Carlo steps that are necessary to 
calculate the relaxation law on the time interval (0, fir n ) is given by 



For the nanoparticle ensembles considered here, Eq. (|3.7J1 for R = 1 and ji = 0.2 yields 
M » 4.33 x 10 10 for T = 300 K, and M « 6.08 x 10 22 for T = 150 K. Such values of M 
render of course the use of that method impractical. For comparison, in our approach the 
number M of time steps At n+1 , defined by the condition J2m=i A£ m = jwr n , equals 157 and 
169, respectively. 

IV. CONCLUSIONS 

We have developed a new method for the numerical simulation of thermally activated 
magnetic relaxation in 2D ensembles of uniaxial ferromagnetic nanoparticles whose easy 
axes of magnetization are perpendicular to their distribution plane. It is based on the 
analytical determination of the probability densities of reorientation of the nanoparticle 
magnetic moments and on the numerical determination of the nanoparticle ensemble states 
for a discrete sequence of times. Using the backward Fokker-Planck equation, we have 
formulated a rigorous approach to calculate those probability densities, and in the case 
of high potential barriers between the equilibrium directions of the nanoparticle magnetic 
moments we have studied the law of magnetic relaxation by this method. 

We have shown that magnetic relaxation in finite nanoparticle ensembles can differ 
strongly from that predicted by the mean-field approximation for infinite ensembles. This 
difference is caused by the finiteness of the ensemble size as well as correlations between 
the magnetic moments, which result from the dipolar interaction between nanoparticles. 
In a finite ensemble, magnetic relaxation for small and large times occurs slower than the 
mean-field theory predicts for infinite ensembles, and for intermediate times the correspond- 
ing relaxation curves, depending on the ensemble size and temperature, can intersect twice. 
Increase of the ensemble size enhances relaxation, and in the limiting case of an infinite 




(3.7) 
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ensemble, magnetic relaxation for small times occurs faster and for large times slower than 
for the mean-field theory. This feature of the relaxation law is caused by the correlation 
effects whose role grows with time. 
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FIG. 1: Schematic representation of the 2D nanoparticle ensemble. 




FIG. 2: Plots of p s im(t) for L = 50 (curve 1), pu m (t) (curve 2), and p m f(t) (curve 3). Inset: The 
same plots and the plot of p S i m (t) for L = 70 (curve 4) for small times. 
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FIG. 3: Plots of xl(*) for L = 50 (curve 1), L = 70 (curve 2), and L = oo (curve 3). Inset: The 
same plots for small times. 
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FIG. 4: Plots of p sim {t) (curve 1) and p mf (t) (curve 2) for L = 50 and T = 150 K. Inset: The 
same plots for small times. 
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